Condensed Matter Physics, 2012, Vol. 15, No 2, 23801: 1-TTo] (3SEBLIIQ9L1E) 
DPI: 10.5488/CMP. 1 5.23801 1 JiJj^XWJBJi 



http://www.icmp.lviv.ua/journal 



On the contact values of the density profiles 
in an electric double layer using density 
functional theoryl* 

LB. BhuiyarP, D. HendersorP, s. SokotowskP 

1 Laboratory of Theoretical Physics, Department of Physics, University of Puerto Rico, 
Box 70377, San Juan, Puerto Rico 00936-8377, USA 

2 Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602-5700, USA 

3 Department for the Modelling of Physico-chemical Processes, Faculty of Chemistry, MCSU, 
20031 Lublin, Poland 

Received September 1 8, 201 1 , in final form October 1 8, 201 1 

A recently proposed local second contact value theorem [Henderson D., Boda D., J. Electroanal. Chem., 2005, 
582, 16] for the charge profile of an electric double layer is used in conjunction with the existing Monte Carlo 
data from the literature to assess the contact behavior of the electrode-ion distributions predicted by the density 
functional theory. The results for the contact values of the co- and counterion distributions and their product 
are obtained for the symmetric valency, restricted primitive model planar double layer for a range of electrolyte 
concentrations and temperatures. Overall, the theoretical results satisfy the second contact value theorem rea- 
sonably well, the agreement with the simulations being semi-quantitative or better. The product of the co- and 
counterion contact values as a function of the electrode surface charge density is qualitative with the simula- 
tions with increasing deviations at higher concentrations. 
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1. Introduction 

One of the more interesting recent developments in the electric double layer research has been the 
advancement of contact value theorems involving the charge profile in a primitive model (PM) planar 
double layer (charged hard spheres moving in a dielectric continuum next to a planar electrode) (see, for 
example, references (lpyi). Such exact conditions, or sum rules as they are often called, are important per 
se in statistical mechanics since they permit unambiguous assessment of various approximate theories 
and hence aid in theoretical development. 

The most famous contact theorem in the double layer literature is the one formulated by Henderson 
and Blum | S], and Henderson, Blum, and Lebowitz (HBL) 1 6] over thirty years ago. It is a condition on the 
contact value of the total density profile in a planar double layer, and for a symmetric valency restricted 
primitive model (RPM) (equisized ions in the PM) planar double layer - the model system of interest in 
this paper, the HBL relation reads 

b 2 

gsumW/2)= lgcoid/2) + g cU {dl2))/2 = a+ — . (1) 



* It is a pleasure to dedicate this paper to Dr. Orest Pizio on the occasion of his 60th Birthday. Douglas Henderson recalls with 
fondness his first meeting with "Don Oresto" in Telavi in the Republic of Georgia in the mid 1980s, where Orest and Myroslav 
Holovko invited him to visit Lviv. During this visit, Orest showed him the city sights, including Stefan Banach's grave and the 
Scottish Cafe, where Banach and his colleagues formulated many important theorems in functional analysis. 
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Here g co , g ctr are the co- and counterion singlet distribution functions, d is the common ionic diameter, 
and a — p/ipk^T) is the bulk osmotic coefficient with p being the bulk pressure, /cb the Boltzmann con- 
stant, and T the absolute temperature. The quantity b = zeo I {eoe x k-&Tx) is a dimensionless parameter 
where z, e, and a are, respectively, the absolute value of the ionic valency, the magnitude of the elemen- 
tary charge, and the uniform surface charge density on the electrode, eo is the vacuum permittivity, and 
e r is the relative permittivity of the continuum solvent. Also, x = \/z 2 e 2 p I {eoe x k^,T) is the Debye-Hiickel 
parameter (inverse Debye screening length) with p-Y.i Pi where p, is the mean number density of the 
ith ionic species. In the RPM case the contact distance, that is, the distance of closest approach of an ion 
to the electrode, occurs at dl2, where d co = d ctI = d. Note that gsumM (x is the perpendicular distance 
from the electrode into the solution) is related to the total density profile 

p{x) = Y j p i {x) = Y j p i g i {x) = pgsumM, (2) 
i i 

where p,-(x) is the singlet density profile of the zth species. It is of interest that the second term in equa- 
tion Q} is just the Maxwell stress. Although equation was obtained from statistical mechanics it is 
consistent with Maxwell's equations. 

Equation (1) is a local expression and the consequent ease of its use has made the HBL contact con- 
dition very appealing in double layer research over the years. For example, the classical Gouy-Chapman- 
Stern (GCS) theory of the double layer satisfies equation but with a = 1, the ideal gas value. Thus, 
for an electrolyte with osmotic coefficient substantially different from unity, the GCS theory can lead to 
appreciable errors especially at low surface charges. 

Sum rules such as equation (1) are useful not only for assessing theories but also for the insight they 
provide. For example, because the coion contact value becomes small at a large surface charge and the 
counterion contact value becomes large, according to this equation, the latter contact value increases as 
square of the surface charge density. On the other hand, the local electroneutrality condition (also a sum- 
rule) requires that the area of the charge profile be equal but opposite in sign to the electrode charge 
density As a result this area increases linearly with the electrode charge density. Consequently, at large 
electrode charge, the oscillations and charge inversions in the charge profile should diminish but cannot 
disappear altogether relative to the contact value. 

Analogous relations for the contact value of the total charge profile in the double layer - the theme 
of the present paper, have been relatively recent. A formal, rigorous relation was derived by Holovko et 
al. dES and Holovko and di Caprio \ 2] using the Bogoliuobov-Born-Green-Yvon hierarchy. For symmet- 
ric valency RPM planar double layer their expression is as follows: 



g diS {d/2) = -zep I dxg sum M— , (3) 

d/2 



where ft = HiksT), y/(x) is the mean electrostatic potential, and gdiffM = [gctrM - gcoM]/2. Again, 
gdiff(x) is now related to the total charge profile, viz., 

q{x) = ej^zipigi (x) = zepg diff (x), (4) 

i 

with Zj being the valency of the ionic species i and z = z co = -z ctr . The definition of gdiffM is, for con- 
venience only, designed to make this quantity positive in general. Since the use of this equation implies a 
knowledge of y/(.x), and giix) throughout the double layer, the expression is non-local. 

Independently, Henderson and Boda (HB) |3J] have proposed an approximate, local expression for 
g(jiff(rf/2) at low electrode charges from empirical considerations, viz., 

g Am {d/2) = ab+0{b 3 ). (5) 

To date a formal, analytic connection between equations l[3) and (5) remains obscure, although in a later 
paper Holovko et al. 1 11] have outlined a very approximate connection. In a series of papers Henderson 
and Bhuiyan |4] and Bhuiyan and co-workers 1 12-15] have tested equation (3) against exact Monte Carlo 
(MC) simulation data for a spectrum of physical states including asymmetric electrolytes fl3Hl5ll and 
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found the equation to be remarkably consistent with the simulations. Theoretical support came from an 
application of the modified Poisson-Boltzmann (MPB) equation, which was found to satisfy equation (3) 
to a very good degree llilllal . Bhuiyan and Henderson (3 have also compared the two relations (equa- 
tions l[2} and l[3)) numerically using simulations and the conclusion was that although exact, equation (2) 
is difficult to implement numerically because of its non-local nature. We note here that an approximate, 
non-local relation for gmf{dl2) has also been suggested by Henderson and Bhuiyan 1 17]. Following con- 
vention, we will call equation {1) the first contact value theorem, while equations (3) and (5) represent 
two versions of the second contact value theorem. 

Another interesting recent result that also concerns us in this paper is the behavior of the product of 
the co- and counterion contact values / = g co {d 1 2) g c x r {d 1 2) in the RPM planar double layer. The classical 
GCS result for this quantity is strictly unity under all circumstances and thus constitutes a basis for the 
classical theory. For example, the value of the counterion contact g ctr (d/2) is as high as the reciprocal 
of the coion contact g co {d/2). However, the corresponding simulation results Qjll dramatically alter the 
classical picture. The product / is seen to be not only different from unity, but also that its characteristics 
as a function of the electrode charge change with the salt concentration. At low concentrations there is a 
maximum before / becomes vanishingly small at high electrode charge. As the concentration increases, 
the height of the maximum decreases and at sufficiently high concentrations the maximum disappears 
completely with / decreasing monotonously. Again, theoretical support for such a behavior of / came 
from the MPB 1 18] and although the hypernetted chain/mean spherical approximation theory does not 
show a maximum, the product / does become very small when the electrode charge is large (l9fl. 

In this study we propose to utilize the HB second contact value theorem and the existing MC simula- 
tion results from the literature for / to assess the density functional theory (DFT) of the planar double 
layer. The DFT has been one of the more successful theories of the electric double layer phenomenon and 
compares favorably with the MPB across planar, cylindrical, and spherical geometries (see for example, 
references 1 20- 22]). Early applications of the DFT to the planar double layer were made by Tang et al. (23| 
and Mier y Teran et al. Hj]. Later Rosenfeld's (Hi techniques were utilized by Mier y Teran et al. ji^] and 
Boda et al. lUH^]. For even recent publications on application of the DFT to the planar double layer, we 
refer the interested reader to the works by Gillespie et al. (HHH], Valisko et al. (3J, Wang et al. [ij. Yu et 
al. (HI], and Pizio et al. liill . Since there is more than one version of the DFT for the planar double layer, 
in the next section we will briefly outline the DFT method used in this paper. Results will be shown in 
section 3, and some conclusions drawn in section 4. 

2. Model and methods 
2.1. Molecular model 

As indicated in the previous section, the model double layer system consists of a binary, symmetric va- 
lency RPM next to a non-penetrable, non-polarizable, uniformly charged planar electrode with a surface 
charge density of o. Since for a given salt concentration, solvent dielectric constant, and temperature, b 
has a linear dependence on a, it is often convenient to specify a in terms of b. 

The ion-ion interaction potential in the Hamiltonian is thus 



where r is the distance between a pair of ions. We also assume that the dielectric constant, e r , is uniform 
throughout the entire system. The bare interaction between an ion of species i and the wall is given by 



where v\ (x) and w (x) are the non-electrostatic and electrostatic (Coulombic) parts of the ion-wall poten- 
tial. The non-electrostatic contribution is a hard-wall potential 




r <d, 
r> d, 



(6) 



Ui{x) = Vi{x)+ Wiix) 



(7) 




x < d/2, 
x > d/2. 



(8) 
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The electrostatic part u>i (x) is given by 



ozie d 

Wi(x) = x, >-. (9) 

e c r 2 



2.2. Density functional theory 



The essence of the density functional theory (see for example, reference |37J]) is that an expression 
for the grand potential, Q, as a functional of the singlet density profiles, pi{x), of each of the species i, 
is initially constructed. At equilibrium the grand potential is minimal with respect to variations in the 
density profiles, viz., 

— = 0. (10) 

opiix) 

This condition is then used to calculate the density profiles and other relevant quantities like the free 
energy. 

In the density functional theory the grand potential of an inhomogeneous fluid can be written in the 
form (iojli^l 

Q = F({p,}) + - Y ezi I pi(x)y(x)dr+ Y / [iii(x) - m]dr, (11) 

/=co,ctr J i'=co,ctr'' 

where p,- denotes the chemical potential of species i. The free energy functions F({p,}) is decomposed into 
ideal (id), hard-sphere (hs), and electrostatic (el) terms as follows F({p,-}) = F-MPi)) + F\ is {{pi}) + F e \{{pi)). 
The ideal term is known exactly 

%«P/})= E fdrtp/Wlnp/W-pit*)]. (12) 

i=co,ctr^ 

For the hard-sphere term, however, we apply the expression resulting from a recent version of the Fun- 
damental Measure Theory 13811 . with the free energy consisting of the terms dependent on scalar and 
vector weighted densities, for details see reference 1391 . 

Following Pizio et al. ^36] electrostatic contribution to the free energy, -F e i({pi}). is represented by 

F e i({p ( }) = f drfMpdx)}), (13) 

where {p; (x)} denotes a set of suitably defined inhomogeneous average densities of a reference fluid. One 
of the simplest possible choices of f e \{{pi (x)}) is to apply the expression resulting from the MSA equation 
of state evaluated via the energy route, namely ji^] 



d ^ 2 _ . . r r 



■A 



f el ({pi(x)})lkT=-— Y z lPi^T-T-j + 



i- co.ctr - L 1 x l ' ^ n 



(14) 



For a symmetric valency situation as in the present case the reduced temperature is T* = 4nk ^J^ e ' d 
Moreover, 

r = (v / l + 2^d-l]/2d. (15) 
The inverse Debye screening length x can be cast in terms of T* 

x 2 = {4nd/T*)Y^Pi{x). (16) 

i 

The last three expressions above correspond to an electroneutral fluid, so that the construction of the 
averaged densities pi{x) at the electroneutrality condition is satisfied. In our approach we follow the 
development proposed by Gillespie et al. j3^|H] described briefly below. 
Let us define the weighted densities p,- (x) as 

PiW = f pi (x')W(lr-r'Ddr', (17) 
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where W(|r- r'|) is a weight function. Gillespie et al. made the assumption, viz., 

WQr-r'\)= ^ — , (18) 

(47r/3)J^(r') 

where 0(|r- r'|) is the step-function. The radius of the sphere over which averaging is performed, Rf, is 
approximated by the "capacitance" radius, that is, by the ion radius plus the screening length 

d 1 

fl f (r)=- + . (19) 

1 2 2T({pi(x)}) 

In addition, Gillespie et al. j3^|H] required that the fluid with the densities {p;(x)} have the same ionic 
strength as the system with weighted densities, {p,(x)}. Consequently, in the case of a symmetric 1:1 
electrolyte the averaged densities {p ; -(x)} are given by 

_ , . _ , . p"i(x) + p 2 (x) 

pi(x) = p 2 (x) = ■ (20) 

Because equations fTs] and (T7J— (20) are coupled, the evaluation of Rf requires an iteration procedure. 
This iteration loop has to be carried out in addition to the main iteration procedure for evaluating the 
density profiles. 

The mean electrostatic potential i[/{x) is determined by the Poisson equation 

^^ = __f_ 2>p ( -(x). (21) 
dx z eoe r i 

The integration of the Poisson equation is carried out subject to the boundary conditions rim z „ 00 i//'(x) = 
andlim^ooi/Ztx) = 0. 

Having specified all the contributions to the free energy functional, the requisite density profiles can 
be obtained by minimizing the grand potential (cf. equation til)). 
All the details of our approach can be found in reference 1 36] . 



3. Results and discussion 

The DFT equations have been solved numerically using the established methods (see for example, 
references |23| 3(J[4(J. We will also present the classical GCS results for comparison purposes, which for 
the RPM case can be obtained analytically. It is convenient to discuss the results in terms of universal 
reduced parameters such as the reduced density p* = pjd 3 and the reduced temperature T* defined 
earlier. Calculations were done at two different reduced temperatures, T* = 0.150 and 0.595, respectively, 
and at each reduced temperature a number of physical states were treated. The value of the ionic diam- 
eter was kept at d = 4.25 xl0~ 10 m throughout. Although a 1:1 valency system was used in the actual 
calculations, in view of universality of T* this becomes a moot point since for a given T* a 1:1 system 
at T is equivalent to a 2:2 system at 4T. For example, in the specific case of T* = 0.15, a 1:1 valency case 
corresponds to -75 K, while a 2:2 valency case corresponds to -300 K. 

In implementing the HB contact condition, Henderson and Bhuiyan |4] found it convenient to recast 
equation (5) in the form 

limf^f^U, (22) 

for symmetrical valency electrolytes. We have followed the procedure here. We note though that a 
straightforward linear plot of equation l[5) with a as the slope has also been done ["16]. In figures [l}[6] 
we present the results for gdm{dl2,b)lb and the contact product / = g co {dl2)g cir {dl2) as functions b 
for T* = 0.15 at p* = 0.02 (c = 0.216 mol/dm 3 ), 0.03 (c = 0.324 mol/dm 3 ), 0.05 (c = 0.541 mol/dm 3 ), 0.10 
(c = 1.08 mol/dm 3 ), 0.20 (c = 2.16 mol/dm 3 ), and 0.25 (c = 2.70 mol/dm 3 ), respectively. The lone filled cir- 
cle on the vertical axis in the upper panel of a figure corresponds to the osmotic coefficient a, which is 
evaluated from the simulations at b = using equation Q). Noticeable immediately from the figures is 
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Figure 1. g di ff(d/2,fc)/Z> (upper panel) and 
g co {dl2)g c xi(dl2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.02 (c = 0.216 mol/dm 3 ) and T* = 
0.15. The symbols represent MC data, while the solid 
line represents the DFT results, and the dash-dotted 
line the GCS results. The filled circle on the vertical 
axis in the upper panel is g sum {d/2,b = 0) = a = 
0.597. MC data from reference 111]. 
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Figure 2. g di ff(d/2,fc)/Z> (upper panel) and 
g co {dl2)g c \ x {dl2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.03 (c = 0.324 mol/dm 3 ) and T* 
= 0.15. The filled circle on the vertical axis in the 
upper panel is g SU m [dl2, b = 0) = a = 0.606. The rest 
of symbols and notation as in figure[l] MC data from 
reference |18||. 
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Figure 3. g^{d/2,b)lb (upper panel) and 
gcoW/2)gctr(d/2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.05 (c = 0.541 mol/dm 3 ) and T* 
= 0.15. The filled circle on the vertical axis in the 
upper panel is g sum [dl2, b = 0) = a = 0.627. The rest 
of symbols and notation as in figure[T] MC data from 
reference fial . 
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Figure 4. g di ff(d/2,fo)/fc (upper panel) and 
gco(rf/2)gctr(rf/2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.10 (c = 1.08 mol/dm 3 ) and T* 
= 0.15. The filled circle on the vertical axis in the 
upper panel is g SU m [dl2, b = 0) = a = 0.684. The rest 
of symbols and notation as in figure[l] MC data from 
reference |'lfj|. 



23801-6 



Contact values of profiles using density functional theory 



the trend that the DFT g^{dl2, b)lb (upper panels of the figures) follows the corresponding simulation 
results very closely for not too high b for the range of concentration treated. Only a very slight discrep- 
ancy is seen at b = 0, which is a consequence of the fact that the DFT does not satisfy the HBL first contact 
theorem exactly. The classical GCS theory satisfies equation l[5) but with a = 1, the ideal gas value. This is 
clear from the figures and for p* = 0.02, 0.03, 0.05, and 0.10 (figures [lj{4), where a is somewhat less than 
unity, the GCS theory leads to deviations from the MC data. The results at a different temperature T* = 
0.595 and at p 3 = 0.00925 (c = 0.1 mol/dm 3 ) and p 3 = 0.0925 (c = 1.0 mol/dm 3 ) are shown in figures[7]and[8] 
respectively. Here too the trends shown by the DFT b)lb and their agreement with the corre- 

sponding simulations are similar to that seen in figures [l]{6] We note that the HNC satisfies equation Q} 
with the first term being a function of the hard sphere compressibility, and equation (5) with a = 1 in the 
first term. For contact values, it is little better than the GCS theory. 

6t 




= 0.25,T =0.15 




Figure 5. g^f{dl2,b)lb (upper panel) and 
g C0 W/2)gctr(d/2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.20 (c = 2.16 mol/dm 3 ) and T* 
= 0.15. The filled circle on the vertical axis in the 
upper panel is g sum (d/2, b = 0) = a = 0.934. The rest 
of symbols and notation as in ngure[TJ MC data from 
reference Hill . 



Figure 6. g^ff{dl2,b)lb (upper panel) and 
gco W/2)g c tr(d/2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.25 (c = 2.70 mol/dm 3 ) and T* 
= 0.15. The filled circle on the vertical axis in the 
upper panel is g sum (d/2, b = 0) = a = 1.09. The rest 
of symbols and notation as in figure[l] MC data from 
reference fl. 



The behavior of the contact product function / is displayed in the lower panel of the figures. Overall 
the characteristics of the DFT plots are in qualitative agreement with the simulations. At p* ^ 0.10, the 
MC data show a maximum with the height of the maximum decreasing as p* increases. The DFT result is 
qualitative and there is a distinct maximum at p* = 0.03 (figure[2}, 0.05 (figure[3), and 0.10 (figure|4), and 
at p* = 0.02, there is the hint of a maximum. Further, the maximum in the DFT curves tends to occur at a 
greater value of b than that in the MC. In figures[5](p* = 0.20) and[7](p* = 0.00925) the plots are initially 
flat, while in figures[6](p* = 0.25) and[8](p* = 0.0925) the initial slope of / is negative. In all of these figures 
the DFT continues to be qualitative with the simulations. The characteristics of the initial slope of / as 
the salt concentration increases can be understood from the following. From equations (l) and (5) one 
has for low b 

g co {d/2)g ctr {d/2) = a 2 + {a-a 2 )b 2 , (23) 

(see for example, equation (18) of reference jl8ll ). This equation is exact in the limit b — «■(). It is easy to 
see at b = that the value of / depends on the value of a. Furthermore, the initial slope of / is negative 
for a < 1 (figures [IHD, the initial slope is approximately zero and the plots are initially flat when a ~ 1 
(figures[5]and[7), and the initial slope is negative when a > 1 (figures[6]and[8). Note again that since in the 
GCS theory a = 1, the right hand side of equation ( fl9t is unity and the initial slope is zero being consistent 
with the observations. 
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Figure 7. g Ai ff(d/2,b)lb (upper panel) and 
gcoW/2)gctr(d/2) (lower panel) as functions of 
6 in a RPM planar double layer for symmetric 
valencies at p* = 0.00925 (c = 0.1 mol/dm 3 ) and T* 
= 0.595. The filled circle on the vertical axis in the 
upper panel is g sum (d/2, b = 0) = a = 0.947. The 
rest of symbols and notation as in Figure 1. MC data 
from reference 11211 . 




Figure 8. g ( ^{dl2,b)lb (upper panel) and 
gco(rf/2)gctr(rf/2) (lower panel) as functions of 
b in a RPM planar double layer for symmetric 
valencies at p* = 0.0925 (c = 1 mol/dm 3 ) and T* 
= 0.595. The filled circle on the vertical axis in the 
upper panel is g sum (d/2, b = 0) = a = 1.31. The rest 
of symbols and notation as in Figure 1. MC data 
from reference 11211 . 



An important property of the simulation data is that the contact product / tends to very small values 
at large values of b. As the surface charge increases the coion population near the electrode is depleted, 
while the counterion population increases. However, the latter also induces packing problems that in- 
hibit distant counterions from migrating too close to the electrode surface. All these lead to the observed 
behavior of /. In the GCS theory however, the decrease in g co is always proportional to the increase in 
g ctr so that classically / = 1 consistently. The DFT / generally follows the MC trend in figures [2j{8] Al- 
though in figure 1 the lack of DFT data beyond b = 8 implies that one cannot be definitive, in view of the 
results in the rest of the figures, it is a fair conjecture that here also the contact product will assume small 
values at still higher values of b. 



4. Conclusions 

In this paper we have examined the predictions of a density functional theory of the planar electric 
double layer with regards to (i) the HB second contact value theorem, and (ii) the behavior of the product 
of the DFT contact values of the co- and counterion distributions vis-a-vis exact MC simulation data from 
the literature. The principal finding regarding (i) is that generally the DFT follows the MC results very 
closely over the range of concentrations and temperatures studied. There is only just a hint of discrepancy 
at b = 0, which is probably tied to the approximation used for the hard-sphere term in the free energy 
functional used to construct the grand potential. By contrast, the GCS results show greater deviations 
from the simulations, especially at lower concentrations when the MC a is less than unity. It is of interest 
to note that the degree to which the DFT satisfies the HB contact condition is very similar to what some of 
us have observed with the MPB theory (l^ll^l with both of the approaches showing slight deviations at b 
= 0. This is not surprising since none of the theories satisfies the HBL first contact value theorem exactly. 

With respect to (ii) above, on the other hand, our calculations indicate that the DFT is broadly in 
qualitative agreement with the characteristics of the simulations, with the product / of the contact values 
tending to small values with increasing surface charge on the electrode. A maximum in / as a function 
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of b seen at lower concentrations although this occurs as, what might be termed, a delayed maximum . 
Importantly though, the behavior of the the initial slope of / as the electrolyte concentration increases 
follows the MC trend. Admittedly, however, there is a quantitative discrepancy between the DFT results 
and the MC data beyond c ~ 1 mol/dm 3 . This is understandable in view of the fact that the product of the 
contact values of the distributions can be a rather more sensitive quantity than their difference so that a 
slight error in either of the contact values tends to become magnified in the contact product. 

The version of density functional that is employed in this paper gives good results for the contact 
values but is less satisfactory in predicting oscillatory profiles. In contrast, other versions of density func- 
tional theory [35,40] are better at producing oscillatory profiles but are less successful for contact values. 
There is more to be done in the development of a fully satisfactory density functional theory. 
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KoHTaKTHi 3HaneHHfl npoc|)miB rycTi/mu b e/iei<Tpi/iHHOMy 
noABiftHOMy tuapi, Bi/iKopi/icTOBytom/i Teopito 
<t>yHKL^ioHa/iy rycTHHH 

/l.B. Eryni-P, fl,. TeHflepcoi-P, c. Coko/iobckP 

/la6opaTopm TeopeTUHHOi cf>i3HKi<i, cf>i3H4HWM cfiaKy^bTeT, YHiBepcuTeT nyepTO PiKO, CLUA 
cPaKy^bTeT xiMfi' i 6ioxiMiT, YHiBepcuTeT BpireMa flHra, npoBO, CLUA 

Biflfli/1 MOfle/llOBaHHfl (fliSWKO-xiMHHUX npOLjeCiB, xiMHHWM c(>aKy/ibTeT, 

YHiBepcuTeT iwi. MapiiCK/iaflOBCbKOi-KK>pi,/lK)6.niH, no/ibma 

HemoflaBHO 3anponoHOBaHy TeopeMy npo /lOKa/ibHe flpyre KOHTaKTHe 3Ha4eHHfl [Henderson D., Boda D., J. 
Electroanal. Chem., 2005, 582, 16] pj\n npocf>i.mo 3ap^fly e/ieKTpn4Horo noflBMHoro luapy noeflHaHO 3 iCHyK)4n- 
mm b /liTepaTypi flaHWMn MoHTe Kap/io 3 MeTOio oi^iHKM KOHTaKTHOi noBefliHKH e/ieicrpofl-iOHHHX po3nofli/iiB, 
nepefl6aneHnx Teopieio cfiyHKujOHa/iy rycTHHH. Pe3y/ibTa™ pj\n koht3kthhx 3Ha4eHb po3nofli/iiB ko- i npoTwio- 
HiB Ta ixHboro flo6yTKy OTpuMaHO fl/ia BwnaflKy cuMeTpuHHOi Ba/ieHTHOcri b paMKax o6Me>KeHOi npuMiTHBHOi 
MOfle/ii n/iocKoro noflBitiHoro luapy fl/w hw3kh KOHLjeHTpaijM i TeMnepaTyp e/ieiapo/iiTy. B uj/iOMy, TeopeTUHHi 
pe3y/ibTa™ .qocnTb flo6pe 3aflOBO/ibHflK)Tb TeopeMy npo flpyre KOHTaKTHe 3Ha4eHHa, y3ro,q>KeHHfl i3 cwMy/ifluj- 
51MI/I - HaniBKi^bKiCHe a6o KpaLue. flo6yTOK ko- i npoTniOHHMX KOHTaKThtux 3Ha4eHb hk cbyHKi4isq rycTWHW 3ap^fly 
noBepxHi e/ieicrpofla AKiCHO y3ro,q>KyeTbCfl 3 cwMy/iflujflMH, a/ie Biflxw/ieHHH Mix o6oMa 3pocrae npw bwluhx 
KOHueHTpaLiinx. 

K/iK>HOBi c/ioBa: eneKTpi/mHui/i noftBtvtHiAvi map, odMexeHa npuMirtiBHa MOftenb, npocpijii rycmHH 
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